arXiv:1503.03096v3 [quant-ph] 11 Mar 2016 


Effective Long-distance Interaction from Short-distance Interaction in Periodically Driven 

One-dimensional Classical System 

Lingzhen Guo 1,2 , Modan Liu 1,3 , and Michael Marthaler 1 
1 Institut fiir Theoretische Festkorperphysik, Karlsruhe Institute of Technology (KIT), D-76131 Karlsruhe, Germany 
2 Department of Microtechnology and Nanoscience (MC2), Chalmers University of Technology, SE-41296 Goteborg, Sweden 
3 Department of Physics, Beijing Normal University, Beijing 100875, China 

(Dated: March 14, 2016) 

We study the classical dynamics of many interacting particles in a periodically driven one-dimensional (ID) 
system. We show that under the rotating wave approximation (RWA), a short-distance ID interaction (5 function 
or hard-core interaction), becomes a long-distance two-dimensional (2D) interaction which only depends on the 
distance in the phase space of the rotating frame. The RWA interaction describes the effect of the interaction on 
the slowly changing amplitude and phase of the oscillating particles, while the fast oscillations take on the role 
of a force carrier, which allows for interaction over much larger effective distances. 
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I. INTRODUCTION 

The properties of one-dimensional (ID) interacting parti¬ 
cles have been investigated as early as 1936 by Tonks 0- 
Later in the 1960’s, Tonk’s model was extended to the more 
general Lieb-Liniger model by allowing the strength of par¬ 
ticles’ contact interaction (6 function interaction) to be arbi¬ 
trary. The exact solution and the thermodynamics of the Lieb- 
Liniger model have been discussed extensively [i2i-[5i]. In the 
last two decades progress in the field of ultracold atoms has 
been signifeant [;6)-l9|] and possibilities to realize a 6 function 
interaction have been proposed EMI- Following these the¬ 
oretical proposals, the Tonks gas and Lieb-Liniger gas were 
observed by several experimental groups lfl3l - [l5n . Due to 
these advances, the properties of interacting particles in one 
dimension have gained renewed interest 11161 1 17TI . Beyond in¬ 
teracting particles in ID free space, the properties of interact¬ 
ing particles in a ID harmonic trap are also of great interest, 
both in theory lfl8l - l25tl and in experiments lfl6L fl7L [26h . 

Ultracold atoms under periodic driving are particularly use¬ 
ful and interesting since periodic driving allows for the gener¬ 
ation of artificial gauge fields li27l - [32ll . Similar to Bloch’s the¬ 
orem for spacial periodicity, the Floquet theory f33l - [37ll can 
be used to treat time-periodic quantum systems. The Floquet 
method transforms a periodic time-dependent Schrodinger eu- 
qation into an eigenvalue problem of a time-independent Flo¬ 
quet Hamiltonian, which is more accessible to a theoretical 
treatment. An intriguing fact of a Floquet Hamiltonian is 
that it can be used to simulate time-independent Hamiltonians 
that are difficult to access otherwise. Additonally the Floquet 
Hamiltonian allows for the study of novel phenomena going 
beyond equilibrium physics 138l - l45ll . In the classical limit, the 
Floquet method corresponds to Poincare mapping I46L147). 
which means we observe the particles every fixed time period 
and mark them on the phase space according to their instan¬ 
taneous positions and momenta. The particles’ trajectories on 
phase space in the discrete time domain can be considered as 
stroboscopic dynamics. Most previous works using Floquet 
theory are based on the single-particle picture, i.e., neglecting 
the interaction between particles. Recent works |4jj|49l] re¬ 



FIG. 1: Sketch of interacting particles in ID harmonic trap: Ul¬ 
tracold atoms (yellow dots) are trapped in ID harmonic potential. 
The atoms are further driven by two counter propagating laser beams 
(red waves). The intensity of laser beams is tuned periodically, i.e., 
oc cos u> d t. The total potential at a fixed moment (blue rippled curve) 
is the harmonic trapping potential plus a cosine function. The inter¬ 
action potential between atoms is assumed to be V(Xj - Xj). 


alized the importance of particles’ interaction in periodically 
driven systems and started to develop scattering theory for 
Floquet states. 

In this work, we investigate the classical dynamics of in¬ 
teracting identical particles in a ID harmonic trap in the pres¬ 
ence of various types of interaction potentials, e.g., 6 func¬ 
tion interaction. Coulomb interaction, hard-core interaction 
and Lennard-Jones interaction. The main goal is to deter¬ 
mine how the particles’ stroboscopic dynamics is influenced 
by their interaction potential V(x, - xj). By going to a rotat¬ 
ing frame and using the rotating wave approximation (RWA), 
we transform the ID spatial interaction V(x, - xj) to a 2D in¬ 
teraction in phase space U(Rjj), which only depends on two 
particles’s distance in phase space Rjj. However, there is a di¬ 
vergence problem in calculating the RWA interaction for, e.g.. 
Coulomb interaction. We analyze the origin of divergence 
and introduce a renormalization procedure to obtain the cor¬ 
rect RWA interaction. Interestingly, we find a short-distance 
(e.g., hard-core interaction) real space interaction can gener¬ 
ate a long-distance RWA interaction on phase space, which 
increases linearly with phase space distance. We justify our 
renormalized RWA interaction by simulating the many-body 
dynamics numerically and compare the results to our theory. 

The article is organized as follows. In section QI] we de¬ 
scribe our model of interacting particles in a ID harmonic trap 



2 


under periodic driving. In section[III] we introduce the RWA 
Hamiltonian by transforming to the rotating frame. We de¬ 
rive the general expression to calculate the RWA interaction 
U(Rjj) for a given real space interaction V(x, - xj). In section 
HY1 we give the canonical equations of motion corresponding 
to both the original Hamiltonian and the RWA Hamiltonian. 
To obtain the stroboscopic dynamics of interacting particles, 
we also dicuss the Poincare mapping method. In section [VI 
we calculate several examples of RWA interactions. We point 
out the divergence problem appearing in the case of Coulomb 
interaction. Then we introduce the renormalization procedure 
to get the correct RWA interaction. We apply our renormal¬ 
ization method to a more general case of an inverse power-law 
interaction potential and the Lennard-Jones interaction poten¬ 
tial. In section[VI] we investigate the two-body dynamics and 
the three-body dynamics with different types of interactions. 
We also justify our RWA interaction by showing the dynamics 
of eight interacting particles under driving. In section [Yll] we 
summarize our work and give an outlook for future work. 


II. MODEL 


The main task of this paper is to identify the role of interaction 
V(x, - Xj) on the slow dynamics of X,(t) and P,(t). 


III. RWA HAMILTONIAN 


We transform the original Hamiltonian 0} into the rotating 
frame with frequency Q./k using the generating function of the 
second kind 
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Here, x = (xi, x 2 , ••■) and P - {P\,P 2 , •■■) represent the assem¬ 
ble canonical coordinates of all the particles. The correspond¬ 
ing canonical transformations of coordinates and momenta are 
given by p, = dG 2 (x, P, t)/dxi, Xj = dG 2 (x, P, t)/dPj, which 
result in the transformation (0. The canonical transformation 
of Hamiltonian itself is given by H(t) + dG 2 /dt. Using the 
rotating wave approximation (RWA), i.e., dropping fast oscil¬ 
lating terms in H(t) + dG 2 /dt, we get the RWA Hamiltonian 
of all interacting particles in the rotating frame. 


We consider many identical particles confined in a ID har¬ 
monic potential and driven by an external driving field as 
sketched in Fig. [j] The particles have the same mass m 
and harmonic frequency u>. For ultracold atoms, the driv¬ 
ing field can be produced bythe interference of two counter- 
propagating laser beams lg,l31 with wavelength a. The inten¬ 
sity of laser beams is tuned periodically, i.e., oc cojoj^t). We 
assume the interaction potential between atoms is V(Xj - xj). 
The total Hamiltonian then is described by 
2 2 

H{t) = Y [ ^- + ^- + A cos(G f)cos(x ; )] + Y V(xi-xj). (1) 

i i<j 

Here, the Hamiltonian has been scaled by the energy unit 
merer/(4n 2 ). The time, coordinate and momentum are scaled 
by or 1 , a/(2n) and nuoaHln) respectively. We also intro¬ 
duced the scaled dimensionless driving frequency D = oj r i/oj 
and driving strength A = An 2 f /{mtJ 2 a 2 ). We define the detun¬ 
ing parameter via 6 = l-Pl/k, where k is a positive integer. In 
this paper, we work in the regime near the resonant condition, 
i.e., [<5| «c 1. Similar Hamiltonians can also be created using 
superconducting devices j50j]. 

We are interested in the case of weak driving regime A <sc 1, 
which means we are not working on a periodic lattice model. 
In our present model, the harmonic trapping potential plays 
an important role and the Hamiltonian 0 does not have spa¬ 
tial periodicity. The basic motion of a particle is dominated 
by the global harmonic oscillation with frequency D./k. The 
driving field and particles’ interactions perturb the phase and 
amplitude of the global harmonic motion. Thus, the total mo¬ 
tion of a particle can be separated into a fast global oscillation 
with frequency Cl/k and a much slower motion representing 
the dynamics of the phase and amplitude of the global motion 

| Xj(t ) = P,{t) sin (ft) + X,(t) cos (ft) 


S = Y | + A cos (y)^(T) cos(fc0,)] + Y ( 3 ) 

i i<j 

Here, we have defined the vector displacement of the z'-th par¬ 
ticle in phase space by = Xj + iP, and two particles’s 
phase space distance by 

R u = \ n e w ‘ - rj e w >\ = ^(Xi-X^ + iPi-Pj) 2 

The detailed derivations are given in the Appendix [0 The 
RWA approximation is valid for small detuning and weak 
driving, i.e., [d[ <s 1 and A <s 1. 

Normally, the original interaction potential V{xi~xj) is only 
a function of the two particles’ distance \xi-xj, which implies 
V(x) - V(-x). Given a particular type of interaction potential 
V(x), we first define the Fourier transformation of the poten¬ 
tial, i.e., V q - jz f j'' dxV(x)e lqx . Then the RWA interaction 
is given by (see more derivation details in Appendix I aT) 

X+OO 

dqVgJoiqRij), (4) 

oo 

where /q(») is the Bessel function of zeroth order. Using the 
integral representation of the Bessel function, i.e., Jo(x) = 
A f ^ e !x sln T dr. we get an alternative form of the RWA in¬ 
teraction as follows 

1 r 2n 2 f? 

U(Rij) = ■ — I V(Rij sin T)dr = — I V(R, j sin t)c/t. (5) 
2?r Jo n J 0 

We see that U(Rjj) is in fact the time average of the interaction 
energy over the oscillation period under RWA. Since the RWA 
interaction U{Rij) is defined in the phase space of the rotating 
frame and it is only a function of phase space distance R,j, we 
call it phase space interaction. 
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IV. CANONICAL EQUATIONS OF MOTION 


The time evolution of the original coordinates, x,(f) and 
Pi(t), of a particle are described by the canonical equations 
of motion (EOM) according to the original Hamiltonian © 

dx, 8H(t) dp, 8H{t ) 

dt dp, ’ dt dxj 

As seen from transformation (0, the values of X,(t) and P,(t) 
can be obtained from the time evolution of x,(t) and p/{t) stro¬ 
boscopic ally every time period of 



(7) 


Here At is defined as the period of the stroboscopic dynamics. 
In this sense, the slow dynamics of Xj(t) and P,(f) is defined 
in the discrete time domain t = mAt with m = 0, ±1, ±2, • • •. 
This technique is called Poincare mapping. 

In the rotating frame, the dynamics of Xj(t) and P,(t) is de¬ 
scribed by the canonical EOM according to the RWA Hamil¬ 
tonian ©, i.e., 




dXj_ = dg_ ^£± = _^S_ 
dt dP, dt dXj 

From the relationship r,e' Si = X, + iPj, we have the explicit 
form of the canonical EOM © as follows (see the details in 
AppendixlBl) 


Xi 


Pi 


6Pi+ 

A cos { k -f )\r k ( ri ) cos^Ofp: - kJ k (n) sin(&0;)f&] 
. ^ ^ P~Pi 


- 6Xj- 

A cos (cos- kJ k (r,) sin(k0,)||] 

-ZjU'(R u )^. 


If we define a complex coordinate via Z, = = A, + iP,, 

the EOM © can be written in an alternative form as follows 


dZj 

dt 


= -[Acos(^) ^~ / ^ r ') sin (k6>,)]z, 


-zjb + A cos (—) k ' cos(k0,)jz, 

2 Vi 

rn Z/ — Z; 


( 10 ) 


We see that the time evolution of Z,-(f) is determined by three 
“forces”. The right-hand-side term in the first line of Eq.dTOli 
is the force parallel to Z, produced by driving. The second line 
of Eg. (Hot is the force perpendicular to Z, produced by detun¬ 
ing and driving. The third line of Eq. dTOl) is the force produced 

z —z 

by the interactions with other particles. The unit vector ^_z \ 
represents the direction from j-th particle to z'-th particle. The 
interaction strength is proportional to the derivative of the 


FIG. 2: Renormalized RWA Interactions: (a) Inverse power-law inter¬ 
action potentials and the example plots for n = 1/2 (Coulomb interaction), 
n = 1, n = 3 and n = oo (Hard-core interaction), (b) Renormalized RWA in¬ 
teraction potentials corresponding to the inverse power-law potentials shown 
in figure (a). Parameter: /( = 0.01. 


RWA interaction, i.e., U'(\Z, - Zj\) = dU(\Zj - Zj\)/d\Zj - Z ; |. 
Different from the interaction in the laboratory, the imaginary 
unit i appearing in this term indicates the direction of this in¬ 
teraction force is perpendicular to the line connecting two par¬ 
ticles. Based on the EOM © or (ITOb . we can calculate all the 
trajectories of interacting particles on phase space. The key of 
the above EOMs is to determine the explicit form of the RWA 
interaction U{Rjj). Below, we will discuss some examples of 
interaction potentials V(r,- rj) and calculate their correspond¬ 
ing RWA interactions U(K,j). 


V. RWA INTERACTIONS 
A. Examples of interaction potentials 

In this subsection, we calculate the RWA interactions for 
several specific interaction potentials U(x, - xj), i.e., the 6 
function interaction potential, the rectangular interaction po¬ 
tential, the hard-core interaction and the Coulomb interaction 
potential. We point out the problem of divergence in the case 
of Coulomb interaction, which can be solved by the renormal¬ 
ization procedure introduced in the next subsection. 


1. 6 function interaction potential 

The 5 function interaction (contact interaction) is used to 
describe the effective interaction between neutral ultra cold 
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atoms in quasi ID confinement liIol - fl5ll . We describe the 6 
function by a Lorenz-function in the limiting case of vanishing 
width, i.e.. 


B e 

V(x, -xj)= /36(xj - xj) = lim -- — - 

£->o n (Xi - Xj)~ + s 2 

We introduce the Lorentz-function here, becuase we will later 
use it for the numerical simulation. Here, fl is the strength of 
the 6 function interaction. The Fourier transformation coeffi¬ 
cient of the above Lorentz function is 14 = fie ^®/(27t). One 
can obtain the RWA interaction from Eq. (0} or Eq. (0 


U(Rii) = - 


4 


]R u + e2 


( 11 ) 


and thus are rebounded back. If the phase space distance of 
two particles is their relative kinetic energy is given by 
EkinfRij) — R~j !4 (see the discussion above Eq. ([27Tb ). For the 
critical condition 17/2 = EkmfRij), Ea. (fl3l) becomes 

U(R,j) ~ -Rtj. ( 14 ) 

n 

As the potential barrier 17 continues to increase beyond the 
critical value /??./2, the RWA interaction (IT4l) should keep un¬ 
changed since the physics in the rest frame does not change 
any more. In the limit of rj — > 00 , the interaction poten¬ 
tial becomes the hard-core potential, i.e., V(x; - xj) = 00 for 
|Xi - JC 7 I < fl and V(xi - Xj) — 0 for |x,- - xj\ > f3. Therefore, the 
formula (IT4l > is the RWA interaction for hard-core interaction. 


In the limit of s —> 0, we have the RWA interaction 


4. Coulomb interaction potential 


U(R,j) = -4-- (12) 

” R ij 

It is interesting to note that the short-distance 5 interaction po¬ 
tential produces an effective long-distance Coulomb-like in¬ 
teraction as function of phase space distance R (/ . 


2. Rectangular interaction potential 

The 6 function interaction is a point-like interaction with 
zero interaction range. Now we allow that the interaction has 
a finite range and define a rectangular interaction potential 

x .- x . f° if \xt - Xj\ >p ; 

77 rect( ' ; ) = l 77/2 if |Xi - Xj\ = [3 ; 

l P [77 if |x,- — Xj\ < p . 

Here, fl is the interaction range and 77 is the interaction 
strength. Applying formula 0. we have (see more details 
in Appendix P 

U(R U ) = /^ arCsin (4) 

\n if Rij < P- 

For R^ » fl, we have the long-range asymptotic behavior 

U(R,j) ~ — 7 J-, for R u »fl. (13) 

7T Rij 

Again, we have an effective long-distance Coulomb-type 
RWA interaction from a short-range real space interaction. 


We now consider the case of Coulomb interaction potential, 
which is the dominant interaction for trapped ions f5ll - (54ll . 
We approach the Coulomb potential by using the following 
function 

V(xj - Xj) = - with s > 0. 

yjiXi - Xj) 2 + S 

Obviously, in the limit of s —> 0, the interaction V(x, - xj) 
becomes the Coulomb potential. Applying formula 0 or 0, 
we obtain (more details provided in AppcndixICl) 


U(R U ) 


2/1 * 1 


71 

Jr 2 . + s 


\ •> 


(15) 


Here, we introduced the effective coupling 


/r s 




R 1 


ij + E 


)■ 


where K(x) is the complete elliptic integral of the first kind. 
Using the asymptotic approximation for the first complete el¬ 
liptic integral [55], i.e., K(ix) ~ - ln(4x) for x » 1 ., Eq. (IT51i 
has the following long-distance asymptotic behavior 

U(Rij) « In (—p), for R u » 1. (16) 

ItKjj v yfs' 

We see that the RWA interaction (IT6l ) diverges in the limit of 
s —> 0, which means the RWA interaction given by Eq. 0 
and Eq.0 is not valid for the Coulomb interaction potential. 
We will analyze the origin of this divergence and introduce 
the renormalization procedure to cancel the divergence. 


3. Hard-core interaction potential 


B. Renormalization procedure 


In the discussion of the rectangular interaction potential, 
we have assumed the potential height 77 is finite. This means 
that the two particles can overcome the potential barrier and 
particles can pass each other if their relative kinetic energy is 
large enough. When 77 is larger than the two particles’ rela¬ 
tive kinetic energy, the particles can not overcome the barrier 


To find the origin of the divergence, we calculate the RWA 
interaction for Coulomb interaction potential using Eq.0 by 
introducing a small cutoff r c as following 


U(R U ) 


= 2^r' 2 ^ dT= W_ ( 2 ), 

TtRij J T _,o Sin T TtRij ''T c / lr c ->0 


(17) 
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We see that U{Rif) diverges in the limit of t c —* 0. Thus, 
the divergence comes from the integral contribution inside the 
small interval [0, r c ]. For a given potential V(r) with r = |x, - 
Xj\, we can estimate the integral inside the interval [0, t c ] by 

2 r Tc 2 r 

U Tc = — I V(R,, sin T)ch ~ - I V(RijT)dT. (18) 
n Jo x Jo 

If U Tc is finite, the potential V(r) is well-behaved inside the 
small distance r < T c Rij. However, for the Coulomb potential 
V(r) oc 1 /r, U T< is divergent. To obtain a finite meaningful 
U ( Rjj ), we subtract this divergence by hand and get the fol¬ 
lowing renormalized RWA interaction 

7T 

2 n 

U(Rjj) = U{Rij) - U Tc = - V(R,j sin j)dr. (19) 
n Jt c 

Equivalently, we we introduce the small cutoff t c to remove 
the divergence. 

The detailed behavior of interaction potential V(r) during 
the collision process is crucial. In the real physical process, if 
V(r) —> oo when r —> 0, the two particles can never touch each 
other since the two particles can not overcome the potential 
barrier. The smallest distance r c they can approach during col¬ 
lision depends on their relative kinetic energy. If the relative 
kinetic energy of two particles is much larger than their inter¬ 
action energy, the collision distance is very short, i.e., r c «: 1 . 
When the two particles are far away, the interaction energy 
can be neglect and their kinetic energy in the center-of-mass 

frame is calculated Efan.(/?,;) = 2 (^) + 4(^r) = Here, 

the kinetic energy of one particle with respect to the center of 
mass is given by ■ We can calculate r c according to the 

energy conservation law in the center-of-mass frame 

V(r c ) = E kin .(R u ) - I Rjj . (20) 

It is important to note here that U(Rjj) is in fact a weighted 
time average of the interaction in the laboratory frame as re¬ 
vealed by Eq. ©. The result is that two particles which col¬ 
lide during an oscillation period, already slow down as the 
they get closer. Therefore, we connect r c and r c by 


Thus, the cutoff r c is given by Eq. (l2Tl) . i.e., t c ~ -p-. There- 

y K ij 

fore, the renormalized RWA interaction for Coulomb interac¬ 
tion is 


2 B 2B* 

U(R,j) = -£-ln(2/r c )= —. (23) 

nRij *Rij 

Here, /J* = f3 ln(/? _I yR4/2) is the renormalized coupling 
strength. The collision factor for Coulomb potential is y — e 2 , 
which is to be calculated in the next section. 


C. Inverse power-law interaction potential 

We now discuss a more general interaction potential form, 
i.e., the inverse power-law interaction potential as plotted in 
Fig. Ha) 


Vto - xj ) 


P ln 

I Xi - Xj\ 2n 


(24) 


We restrict ourselves to integers and half integers n > 1 /2. If 
n — 1 /2, the potential V(x, - xj) has the form of the Coulomb 
potential. If n —» oo, the potential V(x; - Xj) = 0 1 "/\xi - Xj\ 2n 
becomes the hard-core potential with a radius fi. By applying 
Eq. ©, we obtain the RWA interaction 


U(R U ) = 


2p r 12 i 

nRrJ Je—>o sin 2 " r 


W 1 i F( ll 

nR 2 " 1 - 2«l 2 1 z’ 2 



]. (25) 


Due to the term e 1 * 2 ", the above integral diverges in the limit 
of e —> 0 for n > 1/2. Below, we will renormalize U(Rij) for 
integers n > 1, half integers n > 3/2 and n - 1/2 (Coulomb 
potential) respectively. 


1. Integers n > 1 


r c = yRij sin t c ~ yR,jT c (2 1) 

with the important collision factor y > 1 , which is bigger than 
one, because of the extra interaction energy the particles ac¬ 
cumulate as they move closer to each other and start to slow 
down. Therefore we integrate to a radius that can be smaller 
than r c . The collision factor y is phenomenologically intro¬ 
duced here, but we will calculate y for the interaction poten¬ 
tials discussed in this paper. Combining Eq.s ( fl9l >. (l20t and 
(OH . we can obtain the explicit form of renormalized RWA 
interaction U(Rij). 

For example, the collision r c of Coulomb interaction poten¬ 
tial is calculated from Eq. (f20t 


For integers n > 1, we use the properties of special func¬ 
tions 5 - n; | - n; 1 ) = 0 and 2 ^ 1 ( 5 , \ - n\ | - n; 0 ) = 

1. The renormalized RWA interaction can be obtained from 
Eq. (| 25] > by taking e - t c 


U(Rij ) 


2l 2 " 1 l-2n 

7tR 2n 2n - 1 c 

IJ 


(26) 


We use the Eq. < [20b to determine the collision distance r, 


y2n 

' c 


1 ? 

-R 2 : 



(27) 


P_ 

r c 




( 22 ) 


j _j_ 

The truncation r c is given by r c = 2» [iy~ 1 R j . ". Plugging t c 
into Ea. (l26l i. we get the explicit form of renormalized RWA 
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interaction 


U(R U ) 


2/3r"- l 4^' 1-1 
n(2n - 1) 'J 


The collision factor y will be determined later. 


(28) 


2. Half-integers n > 3/2 


For half integers n = k + 1/2 with k > 1, we have the 
following divergence property of Eq. (l25l > 

r ,l 1 3 „ r T(k+\l2)sm{kn + nl2) 

ir\(—,—«;—— --> oo. 

2 2 2 r(k) cos(/br + tt/2) 

The half integers can be approached by taking n — k + 1/2 + e 
with e —» 0. Thus, we have the asymptotic behavior of the 
divergence above. 


Ft 1 1 3 n 

2 ’ 2~ n ’ 2~ 


1 T(k + 1 /2) 
~e V^r(k) 


(29) 


Different from the case of integers n > 1, where we have ne¬ 
glected the zero function iFi{^, \ — n;~ — n; 1) in Ea. (l25l >. the 
divergence (l29b appears in the function 2 ^ 1 ( 5 , \ - n\ f - n; 1) 
for half-integers n > 3/2. It seems that the RWA interaction 
(l28l > is not valid for the case of half integers n > 3/2. How¬ 
ever, we show that this divergence is artificial and is cancelled 
by another divergence in the function of iF\{j, \-n\ |-n;e 2 ). 
To reveal this, we write the function in Taylor’s series (see 
identity (1D4I) in Appendix fDl 


2 ^( 2 . j 


3 2x V 

■n-,--n;e) = 2^ 


m =0 


(l/ 2 ),„(-A-~ c)„, 2m 
*!(!-*-e) m E 


(30) 


The coefficient for m — k in the limit of e —* 0 is 

(1/2 H-k-e) k 1 = 1 F(k + 1/2) 
(1 -k-e\ k\ e sfiY{k) 


This coefficient is divergent as e —* 0 and cancel the diver¬ 
gence of ( l29l >. Therefore, the RWA interaction < 128b is also 
valid for half integers n - k + 1/2 with k > 1. 


3. Coulomb potential n = 1/2 


Assuming n = 1/2 + e/2, the power-law interaction poten¬ 
tial (l24l) becomes V(x,- - xj) = f3\xt - X;|~ 1_e , which goes to the 
Coulomb interaction in the limit of e —> 0. Using the asymp¬ 
totic property, 2 - tv, | —n\ 1 ) » 1 - eln 2 (see the proof 

in Appendix iDl). we have from Eq.(|25]i 



e=0,r c «:l 


(32) 


Compared to formula d23b . it is just the case of Coulomb po¬ 
tential (more details are provided in the Appendix iDb. 


(a) 


(b) 


A 

\ 

\ 

\ 

A 
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-4-20 2 4 
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FIG. 3: Two-body dynamics for different interactions: (a) In the pres¬ 
ence of interaction, the two particles start to rotate at a frequency ojr, which 
depends on their phase space distance R. (b) Rotating frequency cjr as func¬ 
tion of R. The dots represent the data from Poincar^ mapping while the solid 
lines are given by Eq. jh. Different colors represent different interactions as 
indicated on the plot. Interaction potential parameters: ft = 0.1, s = 1 for 
Lorenz interaction; /< 0.1 for Coulomb interaction; e = 0.01. cr = 0.1 for 

Lennard-Jones interaction; Hard-core interaction is modelled by the inverse 
power-law interaction with /? = 0.1 , n = 20. 


4. Collision factor y 

Now we discuss how to determine the collision factor y for 
inverse power-law potential. In principle, the collision fac¬ 
tor y is phenomenologically introduced in Eq.( l2 1 l i. Here, we 
determine it using the correspondence conditions. From the 
expression (l24b we see that the inverse power-law potential 
approaches the hard-core potential in the limit of n —» 00 

P 2n j 00 if \xj - Xj\ < //; 

(xt - Xj ) 2 ' 1 ^ \ 0 if I Xj - xj\ > [3. 

Comparing the RWA interaction ( |28| ) to Eq. d 1 4[ i. we get the 
first correspondence condition 

y 2n ~ 1 !(2n - 1) —> 2 for n —» 00 . (33) 

The simplest assumption is y 2 "~ 1 = 4n + c, where c is a free 
parameter to be further determined. Thus the collision factor 
takes the form of y — {An + c) 1 ^ 2 " -1 ). This form needs to 
be valid for the Coulomb potential, i.e., n - 1/2. By writing 
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n — 1/2 + e/4 with e —* 0, we have 

y — (2 + c + s) 2 ^ e — * (2 + cj t e 2(2+c) for e —-> 0. 

The parameter c can only take the value of -1 to get a mean¬ 
ingful result. Otherwise, the prefactor (2 + c) 2 ^ e —> (2 + c)°° 
takes either zero or infinity, both of which are unphysical. Fi¬ 
nally, we get the expression for the collision factor y for the 
inverse power-law interaction potential 


As see from Eq. (l36l >. the Lennard-Jones potential is com¬ 
posed of two inverse power-law potential with exponents 2m 
and m respectively. During the collision in the range of small 
distance |jc,- - x,j < cr, we have . - . Therefore, 

the term j- in Eq. (l36l > is dominant during the collision. 
The collision factor of the Lennard-Jones potential can be cal¬ 
culated from Ea.(l34l> by choosing n = m, i.e., y = (4m- l)s=r. 


y = (An — l) 2 "->. 


(34) 


VI. MANY-BODY DYNAMICS 


The collision factor becomes y — e 2 for Coulomb interaction 
(n = 1 / 2 ) and y — 1 for hard-core interaction (n —» 00 ). 


A. Two-body dynamics 


5. Summary 


We summarize our results of RWA interaction for the in¬ 
verse power-law potential as following 


U(R U ) = 


4 \n(/3- l yRl/2), 


nR, 

2/3r" 


n(2n- 1) 

%Rn, 


u' 

W.rK 

ij ’ 


for n = j 

for n = 1 , 2 , |, ■ 

for n —» 00 . 


• (35) 


The collision factor is given by y — (4 n - l)s^t. We see that 
the renormalized RWA interaction for the Coulomb potential 
(n = 1/2) still keeps the form of Coulomb’s law, up to log¬ 
arithmic corrections. We show the behaviors of U(R<j) for 
several cases in Fig. [2 b). For every interaction potential with 
n > 1, U(Rij) actually grows with R,j. It is also interesting to 
note that for the case of n — 1, the corresponding C'(R,j) is a 
constant, which means there is no effective interaction in the 
slow dynamics of Xj(t) and Pft). 


6. Lennard-Jones interaction potential 


Another general choice to describe the interaction between 
two noble atoms or molecules is the Lennard-Jones interaction 
potential 


/ cr 2m CT m , 

V(Xi - Xj) = 4e(- --r---—), 

J ' \Xj — X\x: — X:\ m ' 


<J 

I Xj - Xj 


(36) 


where e defines the interaction strength, cr defines the inter¬ 
action range and the parameter m is usually taken rn - 6 in 
the study. By introducing the cutoff r c , we obtain the RWA 
interaction from Ea. (fT9t 



cr 2m 1 
R 2m 2m - 1 


1 — 2 m 
Tr 


cr 


,2m 


_ 1 —m 


R m m - 1 

ij 


)■ (37) 


Then from Ea.d20t and Eq. d2Tli we can calculate the cutoff as 
following 



(38) 


To justify our RWA interaction, we consider the two-body 
dynamics with resonant condition 6 = 0 and zero driving limit 
A -> 0. In Fig. Ha), we show the trajectories of two par¬ 
ticles in phase space with symmetric initial conditions, i.e., 
CXj(O) = \R,Pi(0) = 0) and (X 2 (0) = -\R,P 2 ( 0) = n). With¬ 
out interaction, both particles are performing independent har¬ 
monic oscillations with time period Ar = 27r. Therefore, the 
positions of two particle in phase space are fixed points as 
shown by two black dots in Fig. [3 a). In the presence of in¬ 
teraction, the two particles start to rotate with a frequency o)r 
as shown in the same figure. The rotating frequency u>r de¬ 
pends on phase space distance R of the two particles. In this 
simple case, the motions of the two particles are symmetric. 
The dynamics of either particle, e.g., the first particle, is given 
by Aj(f) = jRcos(u>Rt), P\(t) — jRsin(w«f). The rotating 
frequency cjr can be determined from the EOM @ 


OJR - - 


2 dU(R) 
R dR 


(39) 


Here, U(R ) is the renormalized RWA interaction. 

The two-body rotation frequency cjr can also be obtained 
from Poincare mapping. We simulate the dynamics of two 
particles (xj(f), pi(t)) based on the EOM © and take their val¬ 
ues stroboscopic ally every time period At = 2n. This gives 
the trajectories of two particles in phase space as shown by 
the two dotted lines in Fig. HJa). In Fig. [3jb), we compared 
the numerically extracted ojr to the analytical expression 
(l39b for different types of interactions, i.e., Lorenz interac¬ 
tion, Coulomb interaction, hard-core interaction and Lennard- 
Jones interaction. We see that the result (l39l ) is very good for 
large R and breaks down for small R. From Eg. (PH . we see 
the relative kinetic energy of two particles is R 2 /4. Therefore, 
we conclude that the frequency cor as function of R given by 
EaA3% is valid when the interaction energy is much smaller 
than the relative kinetic energy of two particles, i.e., the cutoff 
given by Eq. (ED satisfies r c <sc 1 . 


B. Three-body dynamics 

As we discussed above, while the two-body dynamics can 
be solved analytically, the three-body dynamics can not be ob¬ 
tained analytically from the RWA EOM (]9} in general. Thus, 
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Coulomb Hard-core Lennard-Jones 



FIG. 4: Three-body dynamics for different interactions: (Left column) The upper figure shows the trajectories of three particles in the presence of Coulomb 
interaction with /? = 0.1. The colored dots represent the initial conditions and the colored solid lines represent the dynamics from RWA EOM The black 

dotted lines are the data from EOM (b) combined with Poincare mapping. The lower figure shows the time evolution of r,-(f) = -yjxf(t) + P 2 (t), i = 1,2,3 
obtained from RWA EOM {9} (solid colored lines) and the Poincare mapping (dotted black lines). (Middle column) The middle two figures show the three-body 
dynamics for hard-core interaction with [3 = 0.1. (Right column) The right two figures show the three-body dynamics for Lennard-Jones interaction with 
parameters e = 0.01 and cr = 0.1. 



FIG. 5: Two-body dynamics under driving: (a) Contour plot of single-particle Hamiltonian function (40) in phase space and the initial conditions of two 
particles (green dots), (b) The trajectories of two particles in phase space without interaction (black, [3 = 0) and with Coulomb interaction (red and blue, = 0.1). 
The black and red trajectories are obtained from the time evolution based on the original EOM {6} combined with the technique of Poincare mapping. The blue 
trajectories are obtained from the time evolution based the RWA EOM 0. Driving parameter: A = 0.1. (c) The trajectories of two particles in phase space with 
interaction potential V(xi — Xj ) = /3 2 /(xi - xj) 2 . All the trajectories with interaction (red and blue, (3 = 0.1 ) and without interaction (black, (3 = 0) completely 
overlap each other. Driving parameter: A = 0.1. 


we solve the three-body problem via numerical simulation. In 
Fig .[4] we compare the trajectories based on the original EOM 
(]6]» with the trajectories given by RWA EOM <|9}. The results 
for the Coulomb, hard-core and Lennard-Jones interactions 
are given in the left, middle and right columns of Fig. [4] re¬ 
spectively. In the upper figure of each column, the colored 
dots (blue, green and red) represent the three particles’ ini¬ 
tial positions in phase space. The colored solid lines are the 
trajectories obtained from RWA EOM (|9). The black dotted 


lines are the results from the original EOM ([6]) combined with 
Poincare mapping. In the lower figure of each column, we 

show the time evolution of r,-(f) = ^jxf(t) + Pj(t), i - 1,2,3 
obtained from RWA EOM (O (solid lines) and the Poincare 
mapping (dotted lines). 

We see that the agreement is good. The discrepancy in the 
long time limit comes from the rotating wave approximation 
we used in this paper. Depending on the type of interactions 
and the initial conditions, the three interacting particles may 
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FIG. 6: Eight-body dynamics under driving: (a) Contour plot of single-particle Hamiltonian function ED in phase space with 6 = 0, A = 0.1 and k = 4. The 
yellow dots indicate the initial conditions of eight particles, (b) The trajectories of eight particles in phase space without interaction (black) and with Coulomb 
interaction (red and blue, = 0.1). Driving parameter: A = 0.1. (c) The trajectories of eight particles in phase space with Lennard-Jones interaction (red and 
blue, cr = 0.1 and e = 0.01). Driving parameter: A = 0.1. 


have complex trajectories in phase space. In the figures for 
Coulomb interaction, particle 2 and particle 3 can be viewed 
as a two-body subsystem while their center of mass forms an¬ 
other larger two-body system with particle 1. The RWA in¬ 
teraction corresponding to Coulomb interaction decays with 
phase space distance. Thus, as shown in the lower figure of 
this column, the rotational frequency of particle 1 around the 
center of particle 2 and 3 is much slower than the rotating 
frequency of particle 2 and 3 around each other. For the fig¬ 
ures for the hard-core interaction, the the corresponding RWA 
interaction increases linearly with the phase space distance. 
Therefore the interaction between particle 1 and particle 2 or 
between particle 1 and particle 3 is stronger than interaction 
between particle 2 and particle 3. As a result, the orbits of 
particle 2 and particle 3 are not as regular as in the case of 
Coulomb interaction due to the strong disturbance by particle 
1. For the Lennard-Jones interaction, we change the initial 
conditions of the three particles. From the lower figure, we 
see that the three particles exchange their positions in phase 
space with a period of about 240At as indicated by the arrows. 
The parameters for the interactions are given in the caption of 

Rg. a 

C. Dynamics under driving 

In the above discussion, we justify the RWA interaction 
without consideration of driving field. Now we add the driving 
term to the EOM and justify the RWA interaction. In Fig. [5] 
we show the two-body dynamics under driving. We choose 
the driving strength A = 0.1 in the RWA Hamiltonian (0. We 
further consider the special resonant condition that the driv¬ 
ing frequency is twice of the harmonic frequency in the RWA 
Hamiltonian 0, i.e., g - g\ + g 2 + U(R\i)- Here, g, with 
i = 1,2 representing the index of two particles is defined as 
the single-particle Hamiltonian function 

gi = + A cos (-n)j k ( ri ) cos(£0,). (40) 


In Fig. 0a) we plot the single-particle Hamiltonian function 
in phase space with chosen parameters <5 = 0 and k = 2, i.e., 
gi = -KJiiri) cos(2 Qj). Due to the driving field, the system 
hosts multiple stable vibrational states with different ampli¬ 
tudes and phases. We only draw the four lowest vibrational 
stable states with different phases (centers of the bright and 
dark regions). Now we put two particles (two green dots) near 
two stable states with phase 0 and n respectively, and study 
their motions in the presence of interaction. In Fig. 0b) and 
Fig. 0c), we compare the full numerical simulations to the 
RWA dynamcis. The trajectories with black and red colors 
are the results from the time evolution based on the original 
EOM 0 with consideration of the two particles’ interaction 
V(xj - Xj) combined with the technique of Poincare mapping. 
The trajectories with blue color are the results from the time 
evolution based on the RWA EOM 0 with the renormalized 
U(Rjj). We see that the results agree with each other. Espe¬ 
cially we find, for the case of n = 1 , the trajectories for /3 = 0 
and f3 ± 0 overlap completely as shown in Fig. 0c). This 
means the inverse power-law potential for n - 1 generates no 
interaction under RWA. 

Fig.0shows the dynamics of eight interacting particles. In 
Fig. 0a) we plot the single-particle Hamiltonian function (l4Qt 
with parameters 6 = 0, A = 0.1 and k — 4. The initial con¬ 
ditions are represented by the eight yellow dots. In Fig. 0b) 
and Fig. 0c) we show the trajectories of eight particles under 
driving. Without interaction the eight particles exhibit local¬ 
ized motions as illustrated by the black curves. In the case 
of Coulomb interaction, the trajectories of eight particles are 
enlarged a bit as shown in Fig. 0b) by the red orbits (from 
Poincare mapping) and blue orbits (from RWA EOM 0). The 
RWA interaction corresponding to Coulomb interactions de¬ 
cays with phase space distance, which can be viewed as small 
perturbation when the particle are far away in phase space. In 
the case of the Lennard-Jones interaction, the trajectories of 
eight particles become global motions as shown in Fig. 0c). 
This reflects the fact that the RWA interaction corresponding 
to the Lennard-Jones interaction increases as the phase space 
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distance increases, which is similar to the hard-core interac¬ 
tion. 


VII. SUMMARY AND OUTLOOK 


We investigated the classical dynamics of periodically 
driven interacting particles in a ID harmonic trap. Under 
RWA, we transform the real space interaction V{x, - xj) into 
a RWA interaction in phase space. Particularly, we find an 
effective long-distance RWA interaction can be produced by 
short-distance real space interactions, e.g., point-like 6 func¬ 
tion interaction and hard-core interaction. The RWA interac¬ 
tion describes the effect of the interaction on the slowly chang¬ 
ing amplitude and phase of the globally oscillating particles, 
while the fast oscillations take on the role of a force carrier, 
which allows for interaction over much larger effective dis¬ 
tances. 

We solved the divergence problem by introducing the renor¬ 
malization procedure. For the Coulomb interaction, our renor¬ 
malization procedure just eliminates the high energy colli¬ 
sion process, which gives rise to a renormalized strength of 
Coulomb interaction. For the hard-core interaction, our renor¬ 
malization procedure gives rise to a completely different long¬ 
distance interaction, which increases linearly with phase space 
distance. We justified our theory by simulating the many-body 
dynamics numerically in the presence of various interaction 
potentials like point-like 6 interaction. Coulomb interaction, 
hard-core interaction and Lennard-Jones interaction. 

The work in this paper only focuses on the classical dy¬ 
namics of interacting ID particles. The next step is to extend 
our study to the quantum regime. An interesting direction is 
combing the lattice structure created by the driving ll23ll . as 
shown in Fig. [6fa), together with the effective RWA interac¬ 
tion. In this way, a Hubbard model can be simulated in phase 
space. This will provide another way to study strongly corre¬ 
lated systems by periodically driving ID system. 


second kind 


G 2 (x,P,t ) = 


XiPi 


x—t cos(Q t/k) 2" 


/ £2 \ 

1 . 

/ O \ 


— p- 

2 ' 

tan (-f). 


(A2) 

we transform to the rotating frame with frequency £l/k. Here, 
x = (xi,x 2 , ■ ■ •) and P = [Pi,P 2 , ■ ■ ■) represent the assemble 
canonical coordinates of all the particles. The corresponding 
canonical transformations of coordinates and momenta is 


Pi = 

which results in 


dG 2 (x., P, t ) 
Ox, 


X, = 


dG 2 (x, P, t ) 
dPi ’ 


(A3) 


P/sin(jr)+A/Cos(yf), p, = P L cos(yt)-A,sin (yt). 


k 
(A4) 


The canonical transformation of H(t ) itself is 


K(t) = H(t) + dG 2 /dt 



+A cos(Or) cos (P t sin (Plt/k) + X ,■ cos(Qf/k)) 

+ ^ V^A Pjj sin(Q t/k) + AX,j cos(Qt/k)) 

«j 



A 
H- 

2 


(cos(Df)e'[ p ' sin < nr A) +x > cos(nr/&)] 



+ ^ V^APjj sin (Ql/k) + AX,j cos (Plt/k)^. (A5) 

i<j 


Here, we have defined the detuning <5 = 1 - Pl/k, the dis¬ 
placement of one particle in phase space r,e' Si = A, + iP,, 
and the relative displacement of two particles in phase space 
(AXjj, APjj) = (Xj - Xj, P, - Pj). Using the Jacobi-Anger ex¬ 
pansion. 
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Appendix A: Derivation of RWA Interaction 

The Hamiltonian of many trapped interacting particles un¬ 
der periodic driving can be written as 


IZCOStp _ 


2 i m J,n{z)e lm \ 


cos(f2r)e' [/>; sin ( nf /U+Y cos(nr/t-)] 
cos(Q t )e ir ‘ co « a,lk - e > ) 

+00 

cos(Dr) i m J m ( ri )e- ,m6 ‘e ,mn,lk 

—00 

1 +00 

~(e in ' + e-‘ n ') i m J m ( ri )e- im0 ‘e imn,/k . 


(A6) 


Hit) = ^ [y + y + AcosfQ t) cos(x,)| + ^ V(x, - Xj). 

i i< i 

(Al) 

Here, all the quantities are scaled to be dimensionless. We are 
working in the regime near the resonant condition, i.e., £2 ~ k 
with a positive integer k. Using the generating function of the 


Under rotating wave approximation(RWA), we drop fast os¬ 
cillating terms in Ea. (lA6l > and obtain 

Cos(Qf)e i J P ; si n(ni/k)+Xi cos(£2//*)] > 

Ui k Jk(rde~ lkBi + r k J-kiri)e Mi ) = e‘ k %J k (n) cos(£0,). 

(A7) 
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Now we focus on the interacting term in Eq. (IA5I) 

V^APjj sin (D.t/k) + A Xjj cos (Clt/kj^ 
= V^Rjj cos (Fit Ik - 9ij)j 

X +oo 

dqV q e iqR u cosmk - ei -' ) 

oo 


+°° /~»+oo 

2 r dqV q J m (qR l j)e ,m(n,lk -^ ) . (A 8 ) 

J-OO 


Here, Rj = (AX,;, A/ J ;/ j, V„ is the coefficient of Fourier 


*■'r ~ ijJ’ r i 

transformation of potential, i.e.. 


v q = 


=i r 

27T J_o, 


dxV(x)e lqx . 


In the RWA, we drop all the terms with m 4 0 in Eq. (IA 8 l > and 
obtain the RWA interaction potential in phase space 

X +oo 

dqV q MqRij). (A9) 

OO 

We then apply the integral representation of Bessel function, 
i.e., 


1 f +lr 

Jo(x) = — J e- xslnT dr, 

and get an alternative form of Eq. (lA9t 

X +OO 

dqV q Jo(qRij) 

OO 

X +oo -I r»+7T 

dqV q — e~ lqR “™ T dT 

OO yj —ft 

i r +n 

= — J V(R,j sm t)c1t. 


(A10) 


Normally, the potential is only a function of the distance [x,- - 
Xj\ which implies V(x) = V(-x). Thus we have 

U(Rij) = ^ J" V(Rjj sin r)dr 

2 fi 

= - I V(Ru sin r)dr. (All) 

* Jo 

Finally, we get the RWA Hamiltonian of all interacting parti¬ 
cles in phase space 


g = ^ [^dr 2 + A cos (y)-4(r ; ) cos(fc0,)j + U(R U ), 
i <<j 


(A12) 


where R// = |r,e ,9 ‘ - r ; e i 6 9| = ^/(X/ - X ; ) 2 + (E, - Py ) 2 is the 
phase space distance. 


Appendix B: Canonical equations of motion in the rotating 
frame 

The canonical equations of motion given by RWA Hamil¬ 
tonian (IA 1 2b in the rotating frame are 


(Bl) 



dXi 

dg 

dPi 

dg 


dt 

dP i ’ 

dt 

dXj' 

From the 

• relationship 



Xi = 

r; cos 6j, 

Pi = n sin 

6, and 

n=xlx 2 + 

we have 






<9(9/ 

-Pi 

ddi 

Xi 


dX t 

x 2 + P y 

dPi 

X 2 + P 2 


dn _ 

Xi 

dn 

Pi 


dXi 


dPi 

yjx 2 + P 2 


(B3) 


As a result, the explicit form of canonical equations of motion 
are 


dX 

dt 


= 5P: 


+A cos ( y)[^(a) cos(n6i)^J- - kJ k (n) sin(fc0,)|^-] 




2 /L *' "dPi 

Pi-Pj 


Pu ’ 


dlA 

dt 


= — SX; 


A cos ( y)[^(a) cos (kOd^- - kJ k (n ) sin(/c0 1 )|^| 




2 ' L ~*' 

Xt-Xj 


Ri 


dXi 

(B4) 


Based on EOM (IB4I) . we can calculate the trajectories of in¬ 
teracting particles in phase space. 


Appendix C: Examples of interaction potentials 
1. Rectangular potential 

We define the following short-range interaction potential by 
a rectangular function 

x .- x . f° if |x, - Xj\ >/9; 

q rect( ; ) = 1 q/2 if |x,- - xj\ = /?; 

Z P 1 // if \xi- Xj \<p. 


Applying formula (IA lib , we have 

2 q r ' 2 ,Rij sin : 


U{R, 


\ 2q j Rij sin r ^ 

reC <—) dT 

Pjj 

4/3?i P v rect(x) 

nR " ^/l - (gx) : 


:dx, (Cl) 
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Rn sinr 


where x = —4^—. Using the definition of rect(x), we have 


Appendix D: Inverse power-law potential 


U(Rij) 


= | T arcsin (|-) if Rij^fr 


if Rij<p. 


For Rjj » // we have the long-range asymptotic behavior 

U(Rij) ~ for R u » (3. (C2) 

n Rn 


2. Coulomb potential 

We approach the Coulomb potential by the following type 
of interaction potential 


V(x, - xj) = 


P 


J(xi - xj ) 2 + e 


with e > 0 . 


The interaction V(x; - xj) goes to Coulomb potential in the 
limit of s —> 0. Applying formula (IA10I) . we obtain 


U(Rij) = 


f 


1 


P 


2n Ir 2 s - n 2 T + E 


2 P 

7i yfs 


f 


w/2 


zdT 


^e _ 1 R 2 . sin 2 r + 1 


-dT 




77 \[s V \[e 

2p* 1 


71 

Jr 2 . + e 


\ V 


(C3) 


Here, we introduced the effective coupling/?* = / 3K [( 

where K(x) is the complete elliptic integral of the first kind 
with the property of 


K(ix) = — =K(J -J--) 

VTT3? v v x 2 + u 


for a real x. We use the asymptotic approximations for the 
first complete elliptic integral 115511 . i.e., K(ix) ~ - ln(4x) for 
x » 1. Then we have the long-range asymptotic behavior of 
potential (1C3I > 


U(R u )*-^-hnC^r), for R u » 1. (C4) 

nRa v V e 


ARi 


We see that in the limit of s —> 0, where the real space inter¬ 
action V(x; - xj) goes to Coulomb potential, the phase space 
interaction potential ( 1C4I ) diverges. We will analyze the physi¬ 
cal origin of this divergence and introduce the renormalization 
procedure to cancel it in the next section 


We assume an interaction potential in the form of inverse 
power-law, i.e., 


V(xj -xj) = 


P : 


> 2/1 


1 |Xj - x ; | 2 " 


(Dl) 


with integers and half integers n > 1 /2. We apply the formula 
(1A lib and obtain 


2 p r 

u <*» - J-l 


2/? r n>1 l 


2 /i 

o sin r 


dr 


2/3 r 1 i 


X 


nRfj Jo t 2n Vl - t 2 

2/? t 


zdt 


3 t x ~ 2n 1 1 3 , it 

t- tF i(—,- n; - n\ r) 

2 " 1 - 2u 2 2 2 2 Uo 


2/3 1 


nR 2 " 


j i r A 1 3 


IJ 

-A-2 n 


2Fl( 2’2 _n; 2" n;£2) LJ- (D2) 

Due to the term e 1-2 ", the above integral diverges in the limit 
of s —> 0 for n > 1/2. We renormalize U(Rjj) for integers 
n > 1, half integers n > 3/2 and n - 1/2 (Coulomb potential) 
respectively. 


1. Identities of hypergeometric function 

To continue, we should introduce some identities of hyper¬ 
geometric function 2 F\{a, b\ c; z) and Gamma function T(z). 

Identity I: Gauss’s theorem|56] of hypergeometric function 

2 F](a,b', c; 1) 

T(c)r(c - a-b) 

2 F\(a,b\ c; 1) = ^-—-p-. /or Re[c] > Re[a + b]. 


F(c - o)F(c - b) 


(D3) 


Identity II: derivatives of hypergeometric function 
2 F\(a, b\ c; z) 

d m , 

- it i (a, b\ c; z) =- 2 r \ (a + m, b + m\ c + nr, z), 

dz m (c) m 


where (#)„, = 


+ m) 

m 


(D4) 


Identity III: Euler’s reflection formula [53] of Gamma func¬ 
tion 


r(l - z)r(z) = — 


. , ■ (D5) 

stn(^z) 

Identity IV: Euler’s duplication formula |57] of Gamma func¬ 
tion 


r(z)r(z + \) = 2 x ~ 2z V^r(2z). 


(D 6 ) 
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Identity V: Taylor’s expansion of Gamma function 


lnT(l+z) = -n + J]^-(-z) k , for |z| < 1, (D7) 


k= 1 


where y ~ 0.5772156649 is the Euler-Mascheroni constant 
and C(k) is the Riemann zeta function at k. Combined with 
Eq.(|D5]». we have the following approximative identity 


T(z) * - e yz , for \z\ 1 , 

z 


2. Half-integers n> 3/2 


(D 8 ) 


For half integers n = k + 1/2 with k > 1, we have the 
following divergence in formula (ID2b 


F(k + 1/2) sin(kn + n/2) 
T(k) cos(Att + 7t/2) 


11 3 

2 Ff-,--n\--n\l) =-yfc 


It seems the phase space potential (128b is not valid for the case 
of half integers n > 3/2. However, we show this divergence 
is artificial and is cancelled by another divergence in formula 
dD 2 l The half integers can be approached by taking n - k + 
1 /2 + e with e —> 0. Then we have 


113 _ 

2^i(-, 2 2 ~” ;1) " A ' 7r 


T(k + 1 /2 + e) cos (kn + Tte) 


T (k + e ) 
1 T(k + 1 12 ) 
e f^F(k) 


sm(kn + ne) 
(D9) 


In fact, there is also a divergent term in the function of 
\ — n\ | - u;s 2 ). To reveal it, we write the function 
iFiO;, 5 - n; | - n; e 2 ) in Taylor’s series using identity (1D4I) 


Ft 1 1 3 

2 ^, 2 ’ 2 - 2 _ £ 


-z 

m —0 
oo 

Z 


2m (l/2) m (-k - e) m r ,l 1 3 

s - 2^1 (-,- n ;- n: 0) 

k\{l-k-e) m 2 2 2 


m =0 


( 1 / 2 )„,(-£-e)„, 2 „, 
k!(l-k-e) m £ 


(DIO) 


We calculate the coefficient for m - k in the limit of e —> 0 
( 1 / 2 \{-k-e) k 1 


(1 -k-e) k k\ 

F(k+l/2) T(-e) F(l-k-e) 


1 


T(l/2) F(-k-e) T(l-e) F{k + 1) 

F(k + 1/2) E(-e) F(l-k-e) 1 

E(l/ 2 ) (-e)r(-e) E(-k-e) r(*+l) 

F(k + 1/2) 1 r(l+k + e) sin(-Gr - 67r) 1 

E(l/ 2 ) ~ T(k + e) sin(;r - kn - en) T(k + 1 ) 

1 E(k + 1/2) 

e V/FT (k) 


(Dll) 


Therefore, this coefficient is divergent as e —> 0 and cancel 
the divergence of ( ID9I ). This means the phase space potential 
is still valid for half integers n = k + 1/2 with k > 1 . 

3. Coulomb potential n = 1/2 

According to identity (ID3b we have 


1 1 3 it 

2Fl( 2’2 _H ’2“"’ 1) 


r (| - n)F{\) 

T(1 -n)T(l) 

T(n) sin(n7r) 


= -yff 


F(n- \) cos {nn)' 


If we assume n - 1 /2 + e/2 with e —> 0, we have 


Ft 1 1 3 

2 -t i (— ,- n: — 

z 1 v 2 2 2 


■ n~, 1 ) 


lF ^2’ 2’ 1 2 ’ 1 "* 
. r(i-f)r(i) 
r(i - f)T(i) 


= yfn 


I’d + §) cos(ne/2) 
T(f) sin(/re/ 2 ) 


1 rq + f) 

f r (f) 

1 r(f)rq + f) 
f r (f ) 2 

2 [ -T(e) 

W 




= 2 ~ 


1 - eln 2 . 


(D12) 


(D13) 


Using the above formula, we have U(Rij) for the potential 
V(xi - xj ) = p\xi - x/l -1-6 , 


26 f- 2n 11 3 , I 1 

U(R U ) = p % n1 0 2 F 1 (-,--n-,--n;t 2 ) . 

/r/? 2 " I -2n 2 2 2 lsinr c 

IJ 

2/3 1 

•u — e mz — r„ ;i 

• c «:l 


(l-6ln2-r- e )| 

+€ -6 \e^,T c < 

i 

2/3 1 


nR l 

u 


, -(1 - eln 2 - 1 + er r lnr c ) 

nR l+e -e U->o,t c «i 

ij 


ln( 2 /r c )| 
nRu U=o, 


i=0,r c «:l 


(D14) 
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